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Abstract 


In some gear dynamic models, the effect of tooth 
flexibility is ignored when the model determines which 
pairs of teeth are in contact. Deflection of loaded teeth 
is not introduced until the equations of motion are 
solved. This means the zone of tooth contact and aver- 
age tooth meshing stiffness are underestimated and the 
individual tooth load is overstated, especially for 
heavily- loaded gears. 

This paper compares the static transmission error 
and dynamic load of heavily-loaded, low-contact-ratio 
spur gears with this effect both neglected and included. 
Neglecting the effect yields an underestimate of reso- 
nance speeds and an overestimate of the dynamic load. 

Nomenclature 


E d 

E. 

E P 

E t 

K « 

Q a .Q b »Q c 

R ai» R a2 

R bl> R b2 

R„ R 2 

S a > S r 


w 


gear error due to tooth deflection, mm 
(in.) 

tooth spacing error, mm (in.) 

tooth profile error or profile modification, 
mm (in.) 

static transmission error of gear pair, mm 
(in.) 

stiffness of gear mesh, N/mm (lb/in.) 

meshing compliance of tooth pair a, b, 
and c, mm/N (in. /lb) 

addendum radii of gear 1 and gear 2, mm 
(in.) 

base radii of gear 1 and gear 2, mm (in.) 

pitch radii of gear 1 and gear 2, mm (in.) 

separation distance in approach and re- 
cess, mm (in.) 

total static transmitted load, N (lb) 


w a ,w b ,w c 

*1 

*2 


static transmitted load on tooth pair a, b, 
and c, N (lb) 

separation angle: rotation of gear 1 

(gear 2 fixed), rad 

separation angle: rotation of gear 2 

(gear 1 fixed), rad 


Subscripts: 

j contact point of meshing tooth pair 

1 driving gear 

2 driven gear 


Introduction 


It is well known that the dynamics of gear systems 
can be influenced considerably by the stiffness of the 
meshing gear teeth. * A principal excitation for gear 
dynamics and vibration is the variation of this stiffness 
caused by teeth entering and leaving mesh. This stiffness 
variation is a primary cause of the time- varying compo- 
nent of static transmission error. The static transmission 
error is defined as relative displacement of the driving 
gear with respect to the driven gear along the line of 
action. The static transmission error can also be affected 
by gear errors such as tooth profile and spacing, runout, 
alignment and deflection under load. 

An important task for developing a gear dynamic 
model is the determination of which pairs of teeth are 
actually in contact at any instant. In some models, the 
gear teeth are treated as rigid when contact conditions 
are determined. 3,4 ' 6 * 9 However, in an actual transmis- 
sion, the load-carrying teeth deform elastically. This 
causes the incoming tooth pair to enter contact earlier 
than the theoretical start of contact. Similarly, the 
loaded outgoing teeth will leave contact later than the 
theoretical end of contact. This extends the tooth con- 
tact zone and increases the average mesh stiffness. 


Copyright • 1993 by H.H. Lin. Published by 
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In this paper, the effect of extended tooth contact 
on heavily-loaded spur gears is examined. The static 
transmission error and dynamic load were calculated for 
gears of moderate contact ratio (1-64) as well as for 
somewhat higher contact ratio (1.95). The calculated 
results were compared to evaluate the influence of 
extended tooth contact on the static and dynamic loads 
of a low-contact-ratio spur gear transmission. The 
findings may form the basis for improvements in the 
spur gear dynamic analysis code DANST. 

Theory and Analysis 

Two sets of low-contact-ratio gears were considered 
for an analytical study. The two sets are the same 
except for the tooth addendum which was adjusted to 
provide contact ratios of 1.64 and 1.95. Parameters for 
the gears are given in Table 1. The analyses were 
performed by the NASA gear dynamic code DANST 
(Dynamic ANalysis of Spur gear Transmissions). The 
analytical procedures are described in the following 
sections. 

Gear System Model 

Figure 1 shows a four-degree-of-freedom, lumped- 
mass model for a typical gear transmission. The model 
includes driving and driven gears, connecting shafts, 
motor, and load. The equations of motion were derived 
from basic gear geometry and elementary vibration 
principles. The dynamic process is studied in the rotat- 
ing plane of the gears and gear tooth contact is assumed 
to be along the line of action. The model and differential 
equations of motion are described in more detail in 
Refs. 10 and II. 

Meshing Stiffness and Transmission Error 
(Neglecting Extended Tooth Contact) 

To study the static transmission error and meshing 
stiffness of a low contact ratio gear system, we designate 
three consecutive tooth pairs a, b, and c, and begin our 
analysis at the moment in which pair a is carrying the 
entire load (single contact zone) and tooth pair b is just 
about to enter contact. The initial contact of tooth pair 
b occurs at the point where the addendum circle of the 
driven gear intersects the line of action. At this instant, 
double contact begins. As the gears rotate, the point of 
contact moves along the line of action. When tooth pair 
b reaches the theoretical point of transition between 
double and single contact, the leading tooth pair a dis- 
engages leaving only pair b in single contact. When 
tooth pair b reaches the next theoretical transition point 
for single and double contact, tooth pair c comes into 
contact and begins to share the load (double contact). 
Thus, the meshing action alternates between double and 
single contact zones as the gears rotate. 


To investigate the effect of tooth flexibility on the 
zone of tooth contact, we will examine in detail the first 
double tooth contact zone (where tooth pairs a and b 
are in contact). With these two tooth pairs in contact, 
the static transmission error E t , and the shared tooth 
load Wj, for each individual tooth pair at contact point 
j may be expressed as: 



w = w* + w b ( 3 ) 

J J 

The tooth spacing errors above are determined with 
reference to tooth pair a (which is therefore assumed to 
have no spacing error). These spacing errors are due to 
manufacturing. The error terms are expressed as linear 
displacements along the line of action. The static trans- 
mission error E t is the total relative displacement of the 
driven gear with respect to the driving gear along this 
line. As long as they are both in contact, the static 
transmission error of tooth pairs a and b must be the 
same. Therefore, from Eqs. (l) to (3), 

q; w ; + (e*), = Q J "w, b + R), + K), < 4 > 

where 

( E .), = M, + M |5) 

( e p), - M), + M, |6) 

M - M), + M, = w > 7 > 

The gear meshing stiffness, K g , at contact point j is then 



In the analyses above and those to follow, the posi- 
tion of the contact point j on the gear teeth is expressed 
in terms of the roll angle of the driving gear tooth. In 
the single contact zone, the transmission error and 
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meshing stiffness equations are much simpler and can be 
derived by similar procedures. 

Gear Teeth Separation Distance 

We define the tooth separation distance as the dis- 
tance between a pair of teeth just out of contact, during 
approach or recess, if there is no elastic deformation. 
This distance, expressed along the line of action, is equal 
to the product of separation angle and base radius of the 
gear. The separation distance will be compared with the 
static transmission error to determine the contact 
condition. 


= cos 


-l 


% 2 ' 




— tan 


-1 


R a2 sin (0 2 + P 2 ) 


(c “ R a2 C0S ^2 ^2) J 

R, 


P 1 




.-1 


'-bi 


= cos 


DO,j 

fRbil 


AO 


(14) 

(15) 

(16) 
(17) 


To calculate the separation distance, we introduce 
the separation angle (exaggerated for clarity in Fig. 2) 
for a pair of teeth (pair b) in approach, where gear 1 
represents the driving gear and gear 2 the driven gear. 
The separation angle is not the same for the two mating 
gears. If gear 1 is held stationary, the separation angle 
A 2 is the angular rotation required for gear 2 to close 
the gap between the teeth of pair b. Likewise, Aj is the 
required rotation of gear 1 while gear 2 is held station- 
ary. The actual tooth contact will start at a point where 
the separation angle of the incoming tooth pair is equal 
to the angular deflection of the preceding tooth pair(s). 

The equations for the separation distance as a func- 
tion of the separation angle can be derived for a tooth 
pair in approach using Fig. 3. In Fig. 3, the driven gear 
(gear 2) is regarded as fixed. Point A represents the 
theoretical start of contact of a tooth pair, point B the 
theoretical end of contact, and point P the pitch point, 
when no load is applied. To find the separation angle 
Aj, the two gears were rotated backward to a position 
just before contact at 6 X and 0 2 , respectively. The tooth 
pair will make contact at point D if the elastic deforma- 
tion of preceding tooth pair(s) caused the driving gear 
(gear 1) to rotate an angle of A r The equivalent separa- 
tion distance along the line of action, Sj, between the 
incoming tooth pair due to the rotation of gear 1 can be 
found from: 


where 


— tan 


^l R bl 

(9) 

Pi ~ “l + *1 

(10) 

— 6 2^2 

(11) 

Pi ~ in v 7l 

(12) 

R»2 sin 02 ] 

1 (is) 


C - R a2 cos p 2 J 


DO, = ^[R a 2 sin (*2 + 0 2 ) f + [ C ~ R a2 cos (®2 + 02 )f 

(18) 


AO, = y[R al sin ^j 2 + [c - R a2 cos ^j 2 


Similar expressions can be derived for the linear separa- 
tion distance S 2 where gear 1 is fixed and gear 2 rotates 
to close the gap. Likewise, expressions can be derived for 
the separation distances Sj and S 2 of a tooth pair in 
recess, where Sj and S 2 are defined for tooth pair a. 

Figure 4 shows the variation of the separation dis- 
tance during approach and recess, as a function of the 
rotation angle (0 X in Fig. 3) for a 1:1 ratio gear pair as 
described in Table 1 and with contact ratio 1.64. The 
zero rotation angle in the abscissa refers to the gear 
position at the theoretical start (or end) of tooth contact 
in approach (or recess). The separation distances Sj and 
S 2 differ and the difference grows larger as the rotation 
angle increases. The magnitude of Sj is less than S 2 in 
approach and greater than S 2 in recess. Since there is no 
particular reason to consider either the driving or the 
driven gear to be fixed, an average value of Sj and S 2 
has been taken as the separation distance. S a is desig- 
nated as the separation distance during approach and S r 
is that for recess. 

The computer program DANST was used to calcu- 
late the static transmission error for the gear system 
described in Table 1. To simplify the analysis, only 
unmodified gears are considered here. DANST is based 
on algorithms developed in Refs. 10 and 11. 

Figure 5(a) shows the theoretical (extended contact 
ignored) static transmission error for gears of contact 
ratio 1.64. The static transmission error is comprised of 
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manufacturing errors added to the deflection due to 
load. In this study, manufacturing errors (such as spac- 
ing and profile errors and runout), are neglected. This is 
a reasonable assumption for high-quality, heavily-loaded 
gears. Therefore, the static transmission error represents 
the elastic deflection of the gear teeth and gear blank. 
The horizontal axis in Fig. 5 is calibrated in terms of 
the roll angle for tooth b. This is the same as the rota- 
tion angle (used in Fig. 4) except for a constant offset. 

Meshing Stiffness and Transmission Error 
(Including Extended Tooth Contact) 



Qj b (s“). + Q‘Q b W 

Of + ? 


( 21 ) 


When E t b , the static transmission error of tooth pair b, 
is less than the separation distance, S r a , tooth pair a 
leaves contact. This occurs as tooth pair b reaches B\ 
(This is the beginning of region III where b is the only 
tooth pair in contact.) The gears remain in single con- 
tact until tooth pair b reaches point C’. 


Superimposed on the transmission error curve in 
Fig. 5(a) are separation distance curves S a and S r . The 
actual point where the approaching tooth pair b makes 
contact is the point labeled C\ where the separation 
distance equals the static transmission error. Likewise, 
tooth pair a, in recess, leaves contact at the point 
labeled B\ 


Five regions (designated I to V) are identified in 
Fig. 5(a). Regions I (AB) and V (CD) represent double 
contact zones; region III (B’C’), is the single contact 
zone; and regions II (BB’) and IV (C’C) represent the 
increased double (or reduced single) contact regions due 
to the effect of tooth flexibility. This effect increases the 
contact ratio of the gear pair about 5 percent (from 1.64 
to 1.72). 


To evaluate the static transmission error of re- 
gion II, we adopt an analysis similar to that presented 
above. We begin at the moment when tooth pair b is in 
contact at the point labeled B in Fig, 5(a). (This is the 
end of the theoretical double contact zone.) Elastic de- 
flection causes tooth pair a to remain in contact until b 
reaches B*. The total transmitted load W shared by 
tooth pairs a and b in this region is: 




( 20 ) 


where the first term at the right hand side of the equa- 
tion represents the load on tooth pair a, and the second 
term represents the load on tooth pair b, at an arbitrary 
contact point j. 

The Btatic transmission error of all tooth pairs in 
contact must be equal, therefore, the transmission error 
of b in this region can be calculated from Ref. 13: 


At the beginning of region IV, when b arrives at 
point C’, tooth pair c engages and gradually increases its 
share of the total transmitted load. The shared tooth 
load and static transmission error of tooth pairs b and 
c change with respect to the rotation of gears. They can 
be determined from the following expressions: 


W = 


_ (<1 ^ M, - fcl 


( 22 ) 


Q" 


Q- 





Q 




(23) 


Figure 5(b) compares the transmission error calcu- 
lated with extended tooth contact included (solid line) 
and ignored (dashed line). 

If the tooth addendum (and hence the height of the 
teeth) is increased, the contact ratio becomes greater. 
The increase in the contact ratio reduces the zone of sin- 
gle contact. The width in the step in the static trans- 
mission error curve will be reduced as the separation 
distance curves S a and S r approach each other. These 
effects can be seen in Fig. 6(a). 

If the theoretical contact ratio is increased to 
slightly less than 2.00, the increase in contact length due 
to extended tooth contact may cause the single contact 
zone to disappear completely. This is illustrated in 
Fig. 6(b) in which the tooth addendum of the gears was 
increased by 20 percent over the standard value to in- 
crease the theoretical contact ratio to 1.95. The actual 
contact ratio (after consideration of tooth flexibility) is 
2.02. 
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Figure 0(a) shows the static transmission error for 
the gears with theoretical contact ratio of 1.95. The 
single tooth contact zone (regions II to IV) is so narrow 
compared to regions I and V that the figure was plotted 
at an expanded scale. (Only portions of regions I and V 
are shown.) For regions I (AB) and V (CD), the static 
transmission error curve is similar to Fig 5(a). Regions 
II (BC’) and IV (B’C) are extended zones of double 
tooth contact (similar to the corresponding regions in 
Fig. 5(a)). The transmission error in these regions can be 
found from Eqs. (18) to (21) as in the previous case. 

In region III (C’B’) tooth pair b is carrying most of 
the load but pairs a and c are also in contact, hence, 
this is a triple contact zone. The transmitted load shared 
by the tooth pairs a, b, and c is given by: 



The transmission error of the three tooth pairs in con- 
tact must be equal, therefore, 

gjfct + folk 

oft* + 

In Fig. 6(b), the magnitude of the transmission 
error was significantly reduced because the single contact 
region was entirely eliminated. The predicted dynamic 
excitation of this gear pair will be similarly reduced 
from that calculated with the extended tooth contact 
neglected. The difference is greater for gears with a 
higher contact ratio (which generally have more flexible 
teeth), especially at heavy load. 

Results and Discussion 



\ * 


(25) 


DANST was used to calculate the dynamic load for 
our sample gear system. To compare the dynamic load 
predicted under different conditions we define a non- 
dimensional term called the dynamic load factor. This 
is the ratio of the maximum dynamic load divided by 
the total static load. The total static load is the torque 
divided by the base circle radius. For gears with contact 
ratio greater than 2, the dynamic load factor may be 
less than 1. 

Figures 7(a) and (b) are comprised of many individ- 
ual solutions for the dynamic load factor arranged in the 
form of speed surveys. The speed surveys are shown with 
the effect of extended tooth contact both neglected and 
included. Figure 7(a) is for a set of sample gears with 
standard tooth addendum (l/DP). The theoretical con- 


tact ratio (neglecting extended tooth contact) is 1.64. 
The response of the gear system peaks at the resonant 
speed near 25 000 rpm. There are also smaller peaks at 
submultiples of the resonant speed. 

Including extended tooth contact in the model 
increases the predicted system resonant speed from 
approximately 28 250 to 24 600 rpm while the predicted 
dynamic load factor at resonance is reduced from about 
2.02 to 1.84, a 9 percent reduction in dynamic load. Ex- 
tended tooth contact results in greater load sharing (in- 
creasing the length of double or triple contact zones) 
which, in-turn, increases the average mesh stiffness. 
Other effects include an increase in the system resonant 
speed and a reduction in the maximum dynamic load. 

Figure 7(b) illustrates the results for a set of 
gears with tooth addendum increased by 20 percent (to 
1.2/DP) over the standard value. This increases the 
theoretical contact ratio to 1.95. Including extended 
tooth contact in the analysis reduces the predicted dy- 
namic factor at resonance from 1.42 to 1.10, a reduction 
of nearly 28 percent. Unlike the example above, there is 
little change in the predicted system resonant speed 
(25 000 rpm). Apparently, elimination of the very nar- 
row single contact zone (Fig. 6) has little effect on the 
average gear meshing stiffness. 

Figures 8(a) and (b) are contour plots which illus- 
trate the effects of both speed and contact ratio on the 
predicted dynamic load factor. The speed was varied 
over the range 2000 to 30 000 rpm and the theoretical 
contact ratio was varied from 1.50 to 1.98. (As above, 
the contact ratio was varied by adjusting the tooth 
addendum.) These figures show how both factors affect 
spur gear dynamics. 

Figure 8(a) shows the predicted results if extended 
contact is neglected. The resonant response at 
23 000 rpm shows the highest dynamic loads for contact 
ratios of 1.52 and 1.70. In Fig. 8(b), the analysis was 
repeated with extended tooth contact included. Fig- 
ure 8(b) shows an overall lower level of dynamic re- 
sponse than Fig. 8(a). The resonant response has shifted 
to about 25 000 rpm, and there is less effect from 
changes in the contact ratio. 

Conclusions 


The NASA gear dynamic code DANST (Dynamic 
ANalysis of Spur gear Transmissions) was used for an 
analytical study of the influence of tooth flexibility to 
extend the zone of tooth contact for heavily-loaded spur 
gears. This effect was both neglected and included as the 
static transmission error and dynamic load was calcu- 
lated for low-contact-ratio spur gears. The following 
conclusions were drawn from the investigation: 
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1. Neglecting the extension of the contact zone 
results in underestimating resonant speeds and overesti- 
mating the dynamic load, especially for heavily-loaded 
gears. 

2. The effect is more significant for gears with a 
theoretical contact ratio nearly (slightly less than) 2.00. 
For these gears, the increased zone of tooth contact may 
extend the actual contact ratio beyond 2.00. 

S. For the cases studied in this paper, ignoring the 
effect results in an underestimate of the contact ratio by 
about 3 to 5 percent. 
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TA BLE 1.— SAMPLE GEAR PARAMETERS 
Number of teeth . . . . * * 


Pressure angle, deg * 

Module, mm (diametral pitch, 1/in.) 3.18 (8) 

Backlash, mm (in.) 0.05 (0.002) 

Face width, mm. (in.) 25.4 (1.00) 

Design torque, N-m (lb-in.) 373 (3290) 

Normalized tooth addendum 1.00 or 1.20 

Theoretical contact ratio * or 1,95 
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(b) Recess. 

Figure 4. — Separation distance calculated for approach and 
recess. 



(a) Theoretical static transmission error and separation 
distance. 



(b) Com pari sk>n of static transmission error with effect 
of extended tooth contact both included and ignored. 

Figure 5.— Static transmission error and separation distance 
for gears with contact ratio = 1 .64. 



.£ 

■o 

c 

.2 


9 > 


(a) Theoretical static transmission error and separation 
distance. (Note: The roll angle axis is plotted at an 
expanded scale so the single tooth contact zone 
may be seen.) 



Extended tooth r 
contact ignored — \ j 


'i 



^Extended tooth 
contact included 


10 1 5 20 25 30 35 

Roll angle, deg 


(b) Comparision of static transmission error with effect 
of extended tooth contact both included and ignored. 


Figure 6. — Static transmission error and separation distance for gears with contact ratio = 1 .95. 
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Dynamic load factor 



(a) Contact ratio = 1 .64. 



Speed, rpm 

(b) Contact ratio = 1 .95. 

Figure 7. — Predicted dynamic load factor of gears as a 
function of speed, with and without the effect of extended 
tooth contact. 



(a) Including effect of extended tooth contact 



2000 11000 20 000 29 000 


Speed, rpm 

(b) Neglecting effect of extended tooth contact. 

Figure 8.— Predicted dynamic load factor of gears as a function of speed 
and contact ratio. 
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